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Abstract. A volume limited and absolute magnitude limited 
sample of A-F type dwarfs within 125 parsecs of the Sun is 
searched for inhomogeneities in the density-velocity space, ex- 
pecting signatures of the cluster evaporation, phase mixing and 
possible disc heating mechanisms. A 3-D wavelet analysis is 
used to extract inhomogeneities, both in the density and ve- 
locity distributions. Thus, a real picture of the phase space is 
produced. Not only are some clusters and streams detected, 
but the fraction of clumped stars can be measured. By esti- 
mating individual stellar ages one can relate the streams and 
clusters to the state of the interstellar medium (ISM) at star 
formation time and provide a quantitative view of cluster evap- 
oration and stream mixing. As a result, we propose a coherent 
interpretation of moving groups or so-called superclusters and 
derive some quantitative evolutionary tracers which we expect 
to serve in the understanding of the large scale evolution of the 
galactic disc. 

The sample is relatively well mixed in the position space since 
less than 7 per cent of the stars are proper motion confirmed 
cluster members. We also detect star evaporation out of the 
Hyades open cluster. 

Only two components of the velocity vectors are provided by 
Hipparcos measurements. Then, the 3D velocity field is re- 
constructed from a statistical convergent point method. The 
wavelet analysis exhibits strong velocity structuring at typ- 
ical scales of velocity dispersion a s tream ~ 6.3, 3.8 and 
2.4 km ■ s _1 . The majority of large scale velocity structures 
(^stream ~ 6.3 km ■ s -1 ) are Eggen's superclusters. As illus- 
trated by the Pleiades supercluster these structures are all char- 
acterized by a large age range which reflects the overall sample 
age distribution. These large velocity dispersion structures rep- 
resent 63% of the sample. This percentage drops to 46% if we 
subtract the velocity background expected by a smooth velocity 
ellipsoid in each structure. Smaller scales (W stream ~ 3.8 and 
2.4 km ■ s _1 ) reveal that superclusters are always substruc- 
tured by 2 or more streams which generally exhibit a coher- 
ent age distribution. At these scales, the contribution of back- 



ground stars is negligible and percentages of stars in streams 
are 38% and 18% respectively. The detailed analysis of the 
phase space structures provides a scenario of kinematical evo- 
lution in the solar neighbourhood: star formation in the galac- 
tic disc occurs in large bursts (possibly subdivided into smaller 
bursts) separated by quiescent periods. The velocity space is 
gradually populated by these star formation bursts which pref- 
erentially fill the center of the velocity ellipsoid. Stars form in 
groups reflecting the clumpy structure of the ISM: about 75% 
of recently formed stars belong to streams whose internal ve- 
locity dispersions do not exceed 4 km ■ s _1 . Most of them dis- 
solve rapidly. A fraction of the initial groups are gravitationally 
bound and form open clusters. Open clusters sustain a longer 
term streaming with quite similar velocity by an evaporation 
process due to internal processes or encounters with permanent 
or transient large mass concentrations. These streams are de- 
tected with a 'stream ~ 2.3 and 3.8 km ■ s _1 and have a coher- 
ent age content. This process explains the survival of streams 
up to 10 9 yr. The existence of streams as old as 2 Gyr seems 
to require other physical mechanisms. The typical scale of so- 
called Eggen's superclusters (Wstream ~ 6.3 km ■ s _1 ) does 
not seem to correspond to any physical entity. The picture they 
form, their frequency and their divisions at smaller scales are 
well compatible with their creation by chance coincidence of 
physically homogeneous smaller scale structures (Jf stream ~ 



3.8 or 2.4 kr 



)• 



Key words: Techniques: image processing, Galaxy: solar 
neighbourhood - open clusters and associations - kinematics 
and dynamics - structure. 



1. Introduction 

1.1. A tentative observation of the stellar gas kinetics 
Providing a complete probe of the kinematics of early type stars 



within a well defined local sphere, Hipparcos data (ESA, 1997) 
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offer the first opportunity to get an unbiased look into the stel- 
lar kinematics at small scales. In paper I ( |Creze et al, 1998 ), we 
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used the average density trend and vertical velocity dispersion 
to map the potential well across the galactic plane. Beyond this 
0-order analysis, we now try to get a picture of the kinematic 
mixture, its inhomogeneities at small scales and their dating. 
This paper explains the methods used and gives the general 
results obtained with our sample. Details both in the methodol- 



ogy and the results can always be found in Chereul et al, 1998 



(hereafter Paper III). A fourth paper will be dedicated to the 
physical interpretation of evolutionary aspects. 
Wavelet analysis is extensively used to extract deviations from 
smooth homogeneity, both in density (clustering) and velocity 
distribution (streaming). A 3-D wavelet analysis tool is first de- 
veloped and calibrated to recognize physical inhomogeneities 
among random fluctuations. It is applied separately in the posi- 
tion space (Section Q) and in the velocity space (Section |]). 
Once significant features are isolated (whether in density or 
in velocity), feature members can be identified and their be- 
haviour can be traced in the complementary space (velocity or 
density). Thus, a real picture of the phase space is produced. 
Not only are some clusters and streams detected, but the frac- 
tion of stars involved in dumpiness can be measured. Then, es- 
timated stellar ages help connecting streams and clusters to the 
state of the ISM at star formation time and providing a quanti- 
tative view of cluster melting and stream mixing at work. 
Only once this picture has been established on a prior-less ba- 
sis we come back to previously known observational facts such 
as clusters and associations, moving groups or so-called super- 
clusters. As a result, we propose a coherent interpretation of 
this phenomenology and derive some quantitative evolutionary 
tracers which we expect to serve in the understanding of the 
large scale evolution of the galactic disc. 
A sine qua non condition for this analysis to make sense is the 
completeness of data within well defined limits. In so far as po- 
sitions, proper motions, distances, magnitudes and colours are 
concerned, the Hipparcos mission did care. Things are unfor- 
tunately not so simple concerning radial velocities and ages. 
More than half radial velocities are missing and the situation 
regarding Stromgren photometry data for age estimation is not 
much better. So we developed palliative methods which are cal- 
ibrated and tested on available true data to circumvent the com- 
pleteness failure. The palliative for radial velocities is based 
on an original combination of the classical convergent point 
method with the wavelet analysis, it is presented in Paper III, 
Section 5.1. The palliative for ages (Section ||) is an empirical 
relationship between age and an absolute-magnitude/colour in- 
dex. It has only statistical significance in a very limited range 
of the HR diagram. 

1.2. Review of known facts 

While clumpy at the time of star formation, the distribution 
of stars in phase space gradually evolves towards a state of 
smoother homogeneity. The precise mechanism of this mix- 
ing is not elucidated: there are reports of moving groups (also 
named superclusters by Eggen) with large velocity dispersions 
(6-8 km ■ s _1 ) surviving after many galactic rotations (Eggen 



1 99 1 , 1 992abc) although not obviously bound by internal grav- 
itation. Considering the very short phase mixing time (a few 
10 s yr, Henon, 1964 1 such streams, exhibiting also large age 
ranges, cannot be explained by the fact that stars originate in a 
same cell of the phase space. 

The study of these moving groups is extensive. A large part of 
the recent work deals with their precise identification among 
early-type (O through F spectral types) main sequence stars 
(Gomez et al, 1990, Chen et al, 1997, Figueras et al, 1997, 
and Chereul et al, 1997) which permit both to probe them 
far from the Sun and to have an individual age estimation 
using Stromgren photometry ( figueras et al, 1991 and Asi- 
ain et al, 1997). Nevertheless, their age content which some- 
times spreads over several hundred million years, as noticed 
in Eggen's work, is still a puzzle. Two main explanations have 
been suggested. On one hand, a supercluster is formed from 
several star formation bursts occurring at different times in 
a common molecular cloud (Weidemann et al, 1992) which 
experiences perturbations such as large scale spiral shocks 
(Comerdn et al, 1997). These bursts can be gravitationally 
bound and maintain streams by cluster evaporation process dur- 
ing their lifetime. Dynamics of star clusters and giant molecular 
clouds (hereafter GMC) in the galactic disc are different. Thus, 
different star formation episodes in a same GMC separated by 
several galactic rotation periods will hardly result in giving a 
supercluster velocity structure. On other hand, asupercluster 
is a chance superposition of several cluster evaporations or 
remnants (and also phase mixing process of unbound recent 
groups) in a same cell of the phase space. Then, a supercluster- 
like velocity structure do not need any physical process to be 
maintained. At a given time in the solar neighbourhood, the 
juxtaposition at random of several cluster remnants creates over 
densities in the velocity field which mimics the existence of a 
physical entity with large velocity dispersion. 



1.3. Hipparcos sample 



Hipparcos data 'ESA, 1997) providing accurate distances and 
proper motions for complete volume limited samples of nearby 
stars, offer the first opportunity to look at the phase mixing 
process and the disc heating mechanisms in action. Such sig- 
natures are searched for with a sample that provides all the 
stars with same physical properties within a well defined vol- 
ume. The sample was pre-selected inside the Hipparcos Input 
Catalogue (ESA, 1992) among the "Survey stars". The limiting 
magnitude is m v < 7.9 + 1.1 • sin | b | for spectral types earlier 
than G5 ( fTuron & Crifo, 1986| ). Spectral types from AO to GO 
with luminosity classes V and VI were kept. Within this pre- 
selection the final choice was based on Hipparcos magnitude 
(m v < 8.0), colours (—0.1 < B — V < 0.6), and parallaxes 
{ir > 8mas). The sample studied (see sample named hl25 in 



Paper I, Creze et al, 1998) is a slice in absolute magnitude of 



this selection containing 2977 A-F type dwarf stars with abso- 
lute magnitudes brighter than 2.5. It is complete within 125 pc 
from the Sun. 
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1.4. Where to find... 

The wavelet analysis procedure is described and discussed 
in Section ^| Age determination methods follow in Section 
[|. Main results of the density analysis in position space 
(clustering) are given in Section [I| The analysis of the velocity 
space (streaming) is given in Section g. A critical review of 
uncompleteness and other systematic effects can be found 
in Paper III, Sections 5.1 and 5.2. Section [T5| presents our 
understanding of the supercluster concept based on all the 
results detailed in Paper III. Conclusions in Section || present 
a simple scenario which organizes the main results of this 
paper to explain observed phase space structures in the solar 
neighbourhood. 



2. Density-velocity analysis using wavelet transform 

An objective method should first be adopted to identify struc- 
tures and determine their characteristic scales and amplitudes. 
The wavelet transform does provide such an accurate local de- 
scription of signal characteristics by projecting it onto a basis 
of compactly supported functions (Daubechies, 1991 ). The ba- 



sis of wavelet functions is obtained by dilatation and translation 
of a unique, oscillating and zero integral function: the mother 
wavelet. The wavelet representation gives signal characteristics 
in terms of location both in position and scale simultaneously. 

2.1. The wavelet transform 

The wavelet transform by ^>(x) of a real one-dimensional 
signal F(x) is defined as the scalar product: 



W s (i) 
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F(x) ■ **( -)dx 



(1) 



where s is the scale and i the position of the analysis. Shape and 
properties of the so-called mother wavelet ^>{x) are similar to 
a Mexican Hat and ensure a quasi-isotropic wavelet transform 
of the signal. It is constructed from a B3 spline function $(x) 
which is compact and regular up to the second order derivative: 

= (| x - 2 | 3 -4 I x - 1 | 3 +6 I a; | 3 



-4 I x + 1 | 3 + I x + 2 | 3 )/12 



(2) 



(3) 



leading to: 

*(.*) = - 

The 3-D scaling function &(x,y,z) used to analyze the ob- 
served F(x,y, z) distribution is a separable function 

$(x,y,z)=$(x)-$(y)-$(z) (4) 

Among several possible implementations of the wavelet analy- 
sis, the "a trou " algorithm, previously used for the analysis of 
large scale distribution of galaxies (Lega et al, 1996, Bijaoui et 
al, 1996) has been selected. The reasons for this choice and the 
principles of the method are given in the following Section. 



2.2. The "a trou" algorithm 

Two main lines of algorithms have been used in order to im- 
plement the wavelet analysis concept. The pyramidal approach 
( Mallat, 1 989| ), used by Meneveau (1991) for turbulence anal- 
ysis, with orthogonal wavelet basis addresses orthogonality 
problems in a rigorous way but it fails providing convenient 
tools to localize structures. On the contrary, the "a trou" 
algorithm associated with non orthogonal wavelet basis, giving 
an identical sampling at each scale (same number of wavelet 
coefficients at each scale), is appropriate for such practical 
purpose. For a detaile d description of the "a trou" algorithm the 
reader is referred to ^olschneider et al (1989) Starck (1993) , 
Bijaoui et al (1996) and details are also provided in Paper III. 
The 3-D distribution F(x, y, z) is binned in a 128 pixel edge 
cube which is wavelet analyzed on five dyadic scales. The 
analyzing function is dilated so that the distance between two 
bins increases by a factor 2 from scale s — 1 to scale s. So 
it is a suitable framework in which local over-densities of 
unknown scales and low amplitude can be pointed at. The 
wavelet coefficient value depends on the gradient of the signal 
F(x, y, z) in the neighbourhood of point (i, j, k) considered at 
a given scale s: the absolute value of the coefficient increases 
all the more since the signal varies on this scale. For a positive 
or null signal like the observed distributions (star counts 
in the solar neighbourhood), over-density structures mainly 
lead to strong positive values of wavelet coefficients while 
under-densities cause strong negative coefficients. 



2.3. Thresholding and segmentation 

Once the wavelet decomposition of the real signal is ob- 
tained, significant wavelet coefficient structures are separated 
from those generated by random fluctuations of a uniform 
background. A thresholding is applied at each scale in the 
space of wavelet coefficients. Thresholds are set at each scale 
and each position by estimating the noise level generated at 
the same scale by a uniform random signal built with the 
same gross-characteristics as the observed one at the position 
considered. Then a segmentation procedure returns pixel by 
pixel the characteristic extent of structures in each dimension. 
The thresholding and segmentation procedures are sketched 
out in Paper III. 



Since we are aiming at an estimation of the total fraction 
of stars actually involved in physical structures, one has to 
care for the casual presence of background stars at the position 
(or velocity) of structures. Cluster membership is classi call y 
tested by the coincidence of tangential velocities (Section 4.1 ). 
This is not possible for stream membership. In this last case, 
only the fraction of non-members can be evaluated (Paper III, 
Section 5.2). Eventually the age distribution of members is 
discussed in the light of various possible scenarii. 
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3. Individual age determination 

3.1. Introduction 

In order to bridge the observed phenomenology with the galaxy 
evolution, it is essential that structures in the phase space be 
dated. Once structure members have been duly identified they 
should be given an age. Stromgren photometry is available for 
some 1608 stars out of a total sample of 2977. When Stromgren 
photometry is available, ages are estimated in a now classical 
way (Section 3.2). There is however a strong suspicion that the 
age distribution of stars observed in Stromgren photometry is 
biased. At least, even though observers are not likely to have 
selected their targets on a prior age indicator, they are likely 
to have favoured stars in clusters and streams, which is highly 
damageable for the present investigation. In order to correct at 
least statistically the suspected resulting bias, we propose an 
empirical palliative age estimation method based on the Hip- 
parcos absolute magnitude and colour. This method is fully de- 
scribed in Paper III. 
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3.2. Ages from Stromgren photometry 

We find 1608 stars with published Stromgren photometry 
(Hauck & Mermilliod, 1990) in our sample. For this sub- 
sample, the effective temperature T e ff, the surface gravity 
log<? and their errors, the metallicity [Fe/H] are derived 
from the photometric indices (b — y), ml, cl, Hp, the 
visual magnitude m v and the rotational velocities v ■ sini 
(if published in Uesugi & Fukuda, 198 1[) us i ng a progr am 



Fig. 1. Stromgren age determinations and their errors for 1077 
stars. 



developed by E. Masana ( Figueras et al, 1991 Masana, 1994 , 
Jordi et al, 1997 ). Based on the above three physical param- 
eters and a model of stellar evolution taking into account 
overshooting effects and mass loss (Schaller et al, 1992), ages 
and masses are inferred with the code developed by Asiain 
(Asiain et al, 1997). By means of the stellar metallicity, this 
algorithm interpolates the set of stellar evolutionary models 
to work at the appropriate metallicity. Ages and masses 
determination were possible for only 1077 stars (a third of 
the sample) due to the failure of the method to get reliable 
metallicities for spectral types between AO and A3. 

The mean error on age determination is 30 % (^0.2 in log- 
arithmic scale) for the bulk of the stars (Figure [l]) and rarely ex- 
ceed 60 % (^0.3 in logarithmic scale). Only the youngest stars, 
between 10 7 and 10 8 years, have errors above 100 % (^0.5 
in logarithmic scale). Nevertheless, this precision is sufficient 
to unambiguously attach them to the youngest age group. The 
distribution of Stromgren ages (Figure ||) for these 1077 stars 
ranges from 10 Myr to 3 Gyr with a peak around 650 Myr 
(log(age) = 8.8). 

3.3. Palliative ages from (M v , (B — V)) 

On a first step, we use existing ages to draw a plot of ages 
versus (Mv, (B-V)). A primary age parameter is assigned as 
the mean age associated to a given range of (Mv, (B-V)). 




Fig. 2. Stromgren age distribution obtained with 1077 stars 
(full line) and palliative age distribution for the same sub- 
sample ( dotted line). 



Then Stromgren age data are used a second time to assign 
a probability distribution of palliative ages as a function of 
the primary age parameter (See details in Paper III). This 
process produces a palliative age distribution which is free 
from possible biases affecting the sub-sample of stars with 
Stromgren photometry for stars with B-V > 0.08. 
The palliative age distribution of stars with Stromgren pho- 
tometry is presented on Figure ^. The small deviations from 
the distribution of original Stromgren ages are due to finite 
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Fig. 3. Stromgren age distribution obtained with 1077 stars 
(full line) and palliative age distribution for the remainder, i.e 
stars without Stromgren photometry (dotted line). 



bin steps used in the discretisation process. The palliative age 
distribution of the sub-sample without Stromgren photometry 
(Figure ||) shows a great difference for very young ages with 
respect to the Stromgren age distribution. The great peak at 
10 7 yr means that more than 400 stars are clearly younger 
than 2-3 10 s yr (see explanations in Paper III) and compose a 
separate population from the rest of the sample. 
These palliative ages have only statistical significance. Hence, 
they sometimes produce artifacts or young ghost peaks in 
some stream age distributions. Such dummy young peaks will 
always appear as the weak counterpart of a heavy peak around 
log(age) — 8.7. Nevertheless palliative ages permit to shed 
light on the age content of the phase space structures when 
Stromgren data are sparse. 



4. Clustering 

4.1. Searching for clusters 

(X,Y,Z) distributions range from -125 pc to +125 pc and are 
binned in a Sun centered orthonormal frame, X-axis towards 
the galactic center, Y-axis in the rotation direction and Z-axis 
towards the north galactic pole. The discrete wavelet analysis 
is performed on five scales: 9.7, 13.6, 21.5, 37.1 and 68.3 pc. 
These values correspond to the size of the dilated filter at each 
scale. 

Over-densities are identified at each scale (Figure by the 
segmentation procedure and the stars belonging to each vol- 
ume are collected. Due to the over-sampling of the signal by 
the "a trou" algorithm, some structures are detected on several 
scales. A cross-correlation has been done between all scales 
to keep the structure at the largest scale provided that there 



is no sub-structure at a lower scale (higher resolution). An 
iterative 2.5 sigma clipping procedure on tangential velocity 
distributions for each group remove field stars and select 
structures with coherent kinematics. 

The same work has been performed to search for under- 
densities (signed by negative wavelet coefficients). No 
convincing evidence of the presence of void has been found. 
In the following we only pay attention to the detection of 
over-densities. 



4.2. Main results 

The space distribution is essentially smooth at all scales. The 
volumes selected by the segmentation procedure contain 10 
per cent of the stars. After the 2.5 sigma clipping procedure 
on the tangential velocities, only 7 per cent are still in clusters 
or groups. Most of them are well known: Hyades, Coma 
Berenices, Ursa Major open clusters (hereafter OC1) and the 
Scorpio-Centaurus association. Otherwise, three new groups, 
probably loose clusters, are detected: Bootes and Pegasus 1 
and 2. Paper III provides a detailed review of all these cluster 
characteristics. Here, we just focus on a newly discovered 
feature concerning the Hyades OC1: its evaporation track. 

The Hyades open cluster's tail is clearly visible at scale 5 
on Figure [I| After selection on tangential velocities, 39 stars 
of the over-density tail are found to have similar motions 
(Figure |]). The tangential velocity component along the 
I axis is nearly the same as the Hyades OCl's one (20.8 
km ■ s" 1 vs. 19.4 km ■ s _1 respectively) but differs along the 
b axis (-2.6 km ■ s" 1 vs. 14.2 km ■ s _1 respectively). The 
extremely peaked age distribution at 5 — 6.3 • 10 8 yr shows that 
stars are slightly younger on average than the Hyades' ones 
(6.3 — 8 • 10 8 yr). Both age determinations are in agreement 
with recent age determination by Perryman et al (1998)| who 



give 6.25 ± 0.5 • 10 yr. Nevertheless, the density distri- 
bution along the tail are a spectacular confirmation of the 
theoretical predictions proposed by Weidemann et al (1992) 
concerning the distributions of evaporated stars from the 
Hyades OC1. These authors expect stars evaporated from the 
Hyades to be distributed in a needlelike ellipsoid centred 
on the cluster center and with longest axis pointing towards 
(I = 97.7°, b = +1.7°). The space distribution of the 39 stars 
exhibit an obvious major axis, with respect to the Hyades 
OC1, pointing towards (7 = 94.0°, b = +7.0°) but are only 
distributed further forward in the direction of Galactic rotation. 
This feature can be produced by escaping stars orbiting closer 
to the Galactic center than the cluster members (Weidemann et 
al, 1992). These stars have a smaller guiding radius implying a 
shorter rotation period. The authors speculate for this type of 
escaping stars, the existence of a phase advance in the vertical 
oscillation which could agree with the difference observed in 
the Tb components. Moreover, we observe these stars higher 
in the plane as it is expected because the open cluster is on the 
upswing. 
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Fig. 4. Significant over-densities in wavelet coefficient spaces for space distributions (grey clumps) and their projections on the 
three main planes (black clumps): (X, Y) plane (bottom) - (Y,Z) plane ( right) - (X,Z) plane (background). Some cluster identifica- 
tions are given. Scale 2, 3, 4 and 5. 



It is striking to notice that stars as massive as 1.8M (which is 
the typical mass in our sample) are also evaporating from the 
cluster. Whatever the cause of the "evaporation" process, slow 
random change of star binding energies due to a large number 
of weak encounters (Chandrasekhar, 1942) or sudden energy 
increase by few close encounters (Henon, 1960) inside the 
cluster, it produces a mass segregation among still clustered 
stars. Indeed, low mass stars are preferentially evaporated and 
massive stars a re retained as s hown by means of nu merical 
simulations in <\arseth (1973) and Ierlevich (1987) and as 
reported in feeid (1992) who find a steeper density gradient 
among the brighter Hyades stars. The asymmetry in the 



density distribution of the escaping stars could have two 
main explanations. We cannot rule out a non detection by 
the wavelet analysis of a symmetric tail provided that its size 
is larger than the coarser scale of analysis. But if it were 
not the case, it might sign an encounter between the Hyades 
OC1 and a massive object on the Galactic center side of 
the cluster. This hypothesis was previously investigated by 



Perryman et al (1998) who find it highly improbable because 
of the velocity with respect to the LSR (30 km ■ s^ 1 ) and the 
mass (~ lO 6 Af ) needed for such an object. A more detailed 
analysis should be done to conclude on the origin of this tail. 
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5. Streaming 

5.1. Mean velocity field 

Global characteristics of the triaxial velocity ellipsoid are 
obtained from a sub-sample of 1 362 stars w hich have observed 
and published radial velocities ( fcSA, 1992 ). The centroid and 
the velocity dispersions in the orthonormal frame centred on 
the Sun's velocity with U-axis towards the galactic center, 
V-axis towards the rotation direction and W-axis towards the 
north galactic pole, are the following: 



- U sample = -10.83 km ■ s _1 and <jxj— 20.26 km ■ s~ 

- V sample = -1 1.17 km ■ s -1 and ay= 12.60 km ■ s~ 

- W sample = -6.94 km ■ s^ 1 and aw= 8.67 km ■ s~ ] 



The sub-sample with observed radial velocities is incomplete 
and contains biases since the stars are not observed at random. 
To keep the benefits of the sample completeness, a statistical 
convergent point method is developed to analyze all the stars 
(see Paper III, Section 5.1). 

5.2. Wavelet analysis of the velocity field 

In a recent paper pehnen (1998) criticizes this methodology 
which could be less rigorous than maximizing the log- 
likelihood of a velocity distribution model. The author argues 
that direct determination of velocity dispersions based on a 
convergent point method produces overestimation and may 
create spurious structures on small scales by noise amplifica- 
tion. Both risks are clearly ruled out by the calibration process 
described above : 
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Fig. 5. The Hyades open cluster's tail after selection on tangen- 
tial velocities (top) and its age distribution (middle), Hyades 
OCl members age distribution (bottom). 



1 . the thresholding is calibrated on numerical experiments so 
as to exclude spurious structure detection at any meaningful 
level (see Paper III, Section 2.2). 

2. the fraction of spurious members in real groups created by 
the convergent point method is estimated from the avail- 
able subset of true 3D velocities based on observed radial 
velocity data (see Paper III, Section 5.2.1). 

3. estimated group velocity dispersions are derived from ra- 
dial velocity subset as well, excluding convergent point re- 
constructed data. 

The confusion in Dehnen's comment is probably related 
to another misunderstanding. In a footnote dedicated to a 
preliminary version of our work, this author explicitly suggests 
that the convergent point method ''which stems for the times 
when more rigorous treatment was impractical on technical 
ground'' was adopted by plain anachronism or ignorance of 
"more rigorous treatment" . Actually, the choice was guided 
by the requirements of the wavelet analysis which gives 
access to a fine perception of 3D structures. Dehnen's analysis 
fits more rigorously a velocity field model with an initial 
coarse resolution of 2-3 km ■ s _1 depending on the velocity 
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Fig. 6. Isosurfaces (grey clumps) in wavelet space after thresh- 
olding for velocity distributions - Projections (black clumps) 
on (U,V) plane (bottom), (V,W) plane (right) and (U,W) plane 
(background) - Scale 1 (3.2 km ■ s^ 1 ). 



component considered while our bin step is 0.8 km ■ s _1 . As a 
result we get (calibrated) significant signal one scale below the 
finest resolution of Dehnen's work. Moreover, the significance 
threshold of velocity field features obtained by Dehnen is 
somewhat arbitrary: features are relevant if they appear in 
more than one of the studied sub-samples. 
Velocities of open cluster stars are poorly reconstructed by 
this method because their members are spatially close. For 
such stars, even a small internal velocity dispersion results a 
poor determination of the convergent point. For this reason, 
we have removed stars belonging to the 6 main identified 
space concentrations (Hyades OC1, Coma Berenices OC1, Ursa 
Major OC1 and Bootes 1, Pegasus 1, Pegasus 2 groups) found 
in the previous spatial analysis. Eventually, the reconstruction 
of the velocity field is performed with 2910 stars. 
Reconstructed (U,V,W) distributions are given in an orthonor- 
mal frame centred in the Sun velocity (see Section 5.1) in 



Fig. 7. Same as Figure ^/or scale 2 (5.5 km, ■ s 1 ). 




Fig. 8. Same as Figure Hfor scale 3 (8.6 km ■ s 1 ). 



the range [-50,50] km ■ s on each component. The wavelet 
analysis is performed on five scales: 3.2, 5.5, 8.6, 14.9 and 27.3 

km ■ s~ x . In the following, the analysis focuses on the first HI, Section 5.2.2 (hereafter procedure 5.2.2). 
three scales revealing the stream-like structures (see Figures 
|(| and ||), larger ones reach the typical size of the velocity 

ellipsoid. Once the segmentation procedure is achieved, stars <- 3 c . , , 

r . j.j. Stream phenomenology 

belonging to velocity clumps are identify. Structures found in 

this reconstructed velocity field are contaminated by spurious 
members created by the method. This contamination is eval- 
uated from the sub-sample with observed Vr by a procedure 
described in Paper III, Section 5.2.1 (hereafter procedure 
5.2.1). The proportion of field stars is also evaluated in Paper 



The A-F type star 3D velocity field turns out highly 
structured at the first three scales. In a previous analysis 
|Chereul et al, 1997 , it has been found that structures are mainly 
revealed in the (U,V) plane rather than in the (U,W) plane. 
This is, probably, the signature of a faster phase mixing pro- 
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cess along the vertical axis. 

In Paper III, Tables 2, 3 and 4 provide mean velocities, ve- 
locity dispersions and numbers of stars remaining after cor- 
rection procedures 5.2.1 and 5.2.2 for streams at respectively 
scale 3, 2 and 1. Here, Table ^summarizes characteristics of the 
streaming organization. For each scale, Table [j] gives the num- 
ber of detected velocity groups after wavelet analysis, the num- 
ber of confirmed streams remaining after elimination of spuri- 
ous members (procedure 5.2.1) and the fraction of confirmed 
stars in streams after corrections 5.2.1 and 5.2.2. The number 
of confirmed streams gives a lower estimation of the real num- 
ber of streams since sometimes the selection of real members 
is performed on very few radial velocity data with respect to 
the potential members. We can notice that the fraction of stars 
involved in streams at scale 2 (~ 38 %) is roughly the same 
as the fraction involved in larger structures at scale 3 (~46%) 
when field stars are removed. It already gives a strong indica- 
tion that large structures at scale 3 could be mainly composed 
by clustering of streams with smaller velocity dispersions from 
scale 2. 

The wavelet scales can be calibrated a posteriori in terms of 
stream velocity dispersions. The distribution of mean velocity 
dispersions (a stream = (j ■ {^u +^ + (7 ™))^) for the streams 
at a given scale yields a better indicator of the characteristic 
scale than the filter size (Figure ^J). The typical velocity disper- 
sions of detected streams at scales 1, 2 and 3 are respectively 
2.4, 3.8 and 6.3 km ■ s -1 . 



Streams appearing at scale 3 (a s t r ea 



6.3 km ■ s 1 ) corre- 



spond to the so-called Eggen superclusters. At smaller scales 
stream ~ 3.8 and 2.4 km ■ s -1 ) superclusters split into dis- 
tinct streams of smaller velocity dispersions. 
The age distribution inside each stream is analyzed. The anal- 
ysis is performed on three different data sets: 



- the whole sample (ages are either Stromgren or palliative), 

- the sample restricted to stars with Stromgren photometry 
data (without selection on radial velocity), 

- the sample restricted to stars with observed (as op- 
posed to reconstructed) radial velocity data (ages are either 
Stromgren or palliative). 



The selection on photometric ages gives a more accurate de- 
scription of the stream age content while the last sample per- 
mits to obtain a reliable kinematic description since stream 
members are selected through procedure 5.2.1. All mean ve- 
locities and velocity dispersions of the streams are calculated 
with the radial velocity data set. Combining results from these 
selected data sets generally brings unambiguous conclusions. 
Results obtained for all the streams found on these 3 scales are 
fully repertoried in Paper III. Here, we focus on the results for 
the Pleiades supercluster (hereafter Pleiades SCI). This exam- 
ple illustrates very well the advance realized in the understand- 
ing of supercluster inner structure. 




Fig. 9. Distribution of stream mean velocity dispersions per 
scale. 



Table 1. Number of streams detected per scale by the wavelet 
analysis and the percentage of stars involved in streams after 
cleaning up spurious members and field stars. 



Scale 


1 


2 


3 


Filter size (km ■ s~ l ) 


3.2 


5.5 


8.6 


stream (klTl • S ) 


2.4 


3.8 


6.3 


Detected velocity groups 
after wavelet analysis 


63 


46 


24 


Confirmed streams after elimination 
of spurious members (Procedure 5.2.1) 


38 


26 


9 


% of stream stars after elimination 
of spurious members (Procedure 5.2.1) 


17.9% 


38.3% 


63.0% 


% of stream stars after elimination 

of background (Procedure 5.2.1 and 5.2.2) 


17.9% 


38.3% 


46.4% 



5.4. Inner structure of the Pleiades supercluster 

The Pleiades supercluster (stream 3-8 in Paper III, Table 2) is 
detected at scale 3 (see Figure [l(] for velocity distributions, Fig- 
1112] for age distributions and Figures |l3| for space dis- 



ures 



tributions) and the set of stars selected on their observed radial 
component gives a mean velocity (U,V,W)=(-14.4,-20. 1,-6.2) 
km ■ and velocity dispersions (cry, cry, cw)=(8.4,5.9,6.3) 
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Fig. 10. Pleiades SCI. Thresholded wavelet coefficient isocon- 
tours at W=-6.2 km ■ s^ 1 of the velocity field at scale 3 
(top), scale 2 (middle) and scale 1 (bottom). Scale 2 shows 
the dichotomy between Cen-Crux and Cen-Lup associations 
(streams 2-26 and 2-12 in Paper III, Table 3) embedded in 
the same clump at scale 3. Scale 1 exhibits clearly the sub- 
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Fig. 11. Pleiades SCI. Age distributions of the Pleiades SCI 
(stream 3-8 in Paper III, Table 2) at scale 3 (top) and scale 
2 (stream 2-5 in Paper III, Table 3) (bottom) . 



km • s _1 . The age distribution (Figure |TT|) covers the whole 
sample age range. The interval is larger than the one mentioned 
by Eggen (1992a): 6-10 6 to 6-10 8 yr. Pure Stromgren age dis- 
tribution, peaks between 6 • 10 8 and 10 9 yr. The Vr selected set 
with palliative ages, shows clearly the preponderance of a very 
young population: a peak at 10 7 yr. A second peak at 6 • 10 s 
yr is also present. Part of this peak can be due to the statistical 
age assignation of very young stars as well as to a real stream. 
But this scale (scale 3) is still too coarse to improve previous 
analysis. 

Scale 2 reveals unambiguously two groups (Figure [IT]) of mean 
ages 10 7 and 10 9 yr in stream 2-5 (Paper III, Table 3). The 
6 • 10 s year old group is no more present among Stromgren age 
distribution and the reminiscence found in the palliative age 
distribution is clearly due to intrinsically very young stars. At 
this scale, the stream is localized at (U,V,W)=(-12.0,-21.6,-5.3) 
km ■ s _1 with velocity dispersions (cry, cry, cw)=(5.3,4.7,5.9) 
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Fig. 12. Pleiades SCI. Age distributions of the two sub-streams: 
stream 1-6 and 1-7 (in Paper III, Table 4) discovered in the 
Pleiades SCI at scale 1. The oldest population is mainly present 
in the stream l-6( top ) while the youngest one composes stream 
7-7 (bottom). 
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Fig. 13. Space distribution of Pleiades SCI obtained from the 
Vr selected sub-sample at scale 3 (stream 3-8 in Paper III, 
Table 3). 



The relative large velocity difference with the Pleiades OC1 (~ 
8 km ■ s _1 ) and the space concentration (Figure |l4|) suggest 
that these stars were formed from the same interstellar cloud 
complex and at the same time as the Pleiades OC1, yet sepa- 
rately. The second stream is much older (roughly 10 9 yr) and 
is probably related to an old cluster loosely bound by internal 
gravitation which is finishing to dissolve now. The probability 
to find such a coincidence in the velocity volume cov ered by 
the Eggen's supercluster is quite high (see Section 5.5 ). 



km ■ s _1 . 

At the highest resolution (scale 1) the main velocity clump 
splits into two components. The first one, at (U,V,W)=(-13.1,- 
21.9,-7.1) km ■ with (ajj, ay, <rw)=(3. 1,3.3,2.5) km ■ s _1 
(stream 1-6 in Table 4, Paper III) contains almost all the old- 
est stars (10 9 yr). The second component at (U,V,W)=(-1 1.1,- 
21.9,-5.9) km ■ s" 1 with (a v , oy, aw>( 1.7,3 .0,1. 9) km ■ s" 1 
(stream 1-7 in Paper EE, Table 4) is almost exclusively com- 
posed of the youngest population (Figure [T^). Moreover, the 
youngest component has significantly smaller velocity disper- 
sions and its space distribution is more clumpy (Figure |l4| ). 
The interpretation comes out quite naturally. The Pleiades SCI 
is a chance superposition of two main streams originating from 
two different star formation epochs. One very young stream is 
a few 10 7 year old and is related with the Pleiades open cluster. 



5.5. How about Eggen 's superclusters 

The review of age and space distributions of the largest veloc- 
ity structures detailed in paper III puts some classification in the 
understanding of the so-called superclusters. As it is shown in 
Section 5.4 for the Pleiades SCI, it is clear that most superclus- 
ters, when looked at sufficient resolutions (small scales), split 
into sub-components covering much restricted age ranges. The 
question then arises whether superclusters are real significant 
structures involving stars born at different epochs but tied to 
the same cell of the phase space by any binding mechanism or 
are they only chance coincidences between streams at smaller 
scales, each essentially coeval. In the first case, a complex sce- 
nario such as the one proposed by Weidemann et al (1992) 
might be necessary. In such a scenario, cluster stars are formed 
at different epochs out of a single molecular cloud and remain 
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Fig. 14. Space distribution of the 2 sub-streams of Pleiades 

SCI from the Vr selected sub-sample at scale 1: stream 1-6 
(top,) and stream 1-7 (bottom) in Paper III, Table 4. 



gravitationally bound to the cloud. A process which looks hard 
to maintain over long periods. 

It is easy to evaluate the probability that superclusters are 
formed by chance coincidence of smaller velocity dispersion 



Table 2. Observed number of superclusters including n or more 
elementary streams, compared with what would be expected if 
38 independent elementary streams were distributed randomly 
in the whole velocity volume. 



n 


Probability 


Probability 


Expected number 


Observed number 








of superclusters 


of superclusters 




k=n 


k>n 


k>n 


k>n 





0.296 


1 






1 


0.360 


0.704 


22.0 


26 


2 


0.219 


0.343 


10.7 


6 


3 


0.089 


0.124 


3.9 


3 


4 


0.027 


0.035 


1.1 


1 


5 


0.007 


0.008 


0.3 


1 


6 


0.001 


0.002 


0.1 






streams. In Paper III, Section 5.2.2 we computed the fraction 
of a smooth gaussian 3D distribution embedded in the veloc- 
ity volume covered by superclusters to estimate the fraction of 
field stars. This fraction was found to be 0.192 for the 6 identi- 
fied superclusters. That is an average of 0.032 per supercluster. 
Under the assumption that the 38 streams appearing at scale 1 
with a stream ~2.4 km ■ s _1 (see Table |l| and Paper III, Ta- 
ble 4) correspond roughly to 38 real independent causes, most 
originating at different epochs and different places, the aver- 
age number of such low dispersion streams expected to occur 
in a typical supercluster velocity volume is 38-0.032=1.216. In 
Table ^ we give the poissonian probability to get 0, 1,2... coin- 
cidences in any supercluster velocity volume, and the number 
of superclusters built out of 1, 2,... elementary streams that one 
would expect to get by chance in the explored velocity vol- 
ume. The total velocity volume contains 1/0.032=31.2 typical 
supercluster volumes. These average numbers should be com- 
pared with the observed statistics of supercluster richness in 
elementary streams. Clearly enough no meaning full test at any 
reasonable significance level can be built to reject the null hy- 
pothesis that superclusters have no physical reality. 
So superclusters most likely result from the chance coinci- 
dence in a large cell of the velocity space of several small 
streams, and the physical interpretation has to be searched for 
only at the smaller scale. Almost all the phenomena observed 
here can be explained by a single scenario resulting from two 
dominant mechanisms: phase mixing and cluster evaporation 
or disruption. At formation most stars form in clump of the 
ISM generating short lived streams which dissolve essentially 
over mixing time scales (~ 10 8 yr): it is probably the case of the 
Centaurus associations and the very young component of the 
Pleiades SCI. Only streams massive enough to create some self 
gravitationally bound systems (more or less loose clusters) sur- 
vive and create, as they dissolve, moderately old streams with 
age between 5-10 8 and 10 9 yr (see streams in Hyades SCI, Sir- 
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ius SCI and New SCI in Paper III). In some cases those mod- 
erately old streams can be explicitly connected to the cluster 
they are escaping from (Ursa Major OC1 for the old compo- 
nent of Sirius SCI and possibly Coma Berenices OC1 for the 
New SCI 610 8 year old component). However, there are a few 
much older groups unrelated to heavy superclusters (streams 
3-9, 3-11 and 3-14 in Paper III, Table 2). In those rare cases, 
all connected to highly eccentric orbits, a completely different 
mechanism should be advocated which might be trapping on 
resonant orbits generated by the potential of the bar (Dehnen, 
1998). 

6. Conclusions 

A systematic multi-scale analysis of both the space and 
velocity distributions of a thin disc young star sample has been 
performed. 

The sample is well mixed in position space since no more 
than 7% of the stars are in concentrated clumps with coherent 
tangential velocities. In this paper we focus on the detection 
of the evaporation of relatively massive stars (1.8 MQ) out of 
the Hyades open cluster. The mapping we realized show an 
asymmetric pattern further forward the Hyades OC1 orbit. This 
could be the signature of a violent encounter on one side of the 
open cluster with a massive molecular cloud. Such a picture 
is in agreement with the origin of streams from open cluster 
evaporation or disruption. 

The 3D velocity field reconstructed from a statistical con- 
vergent point method exhibits a strong structuring at typical 
scales of a 'stream ~ 6.3, 3.8 and 2.4 km ■ s _1 . At large scale 
(scale 3) the majority of structures are identified with Eggen's 
superclusters. These large scale velocity structures are all 
characterized by a large age range which reflects the overall 
sample age distribution. Moreover, few old streams of ~ 2 
Gyr are also extracted at this scale with high U components 
towards the Galactic center (see Paper III). Taking into account 
the fraction of spurious members, evaluated with an observed 
radial velocity data set, into all these large velocity dispersion 
structures we show that they represent 63% of the sample. 
This percentage drops to 46% if we remove the velocity 
background created by a smooth velocity ellipsoid in each 
structure. Smaller scales (cfstream ~ 3.8 and 2.4 km ■ s _1 ) 
reveal that superclusters are always substructured by 2 or more 
streams which generally exhibits a coherent age distribution. 
The older the stream, the more difficult the age segregation 
between close velocity clumps inside the supercluster velocity 
volume. At scale 2 and 1, background stars are negligible and 
percentages of stars in streams, after evaluating the fraction of 
spurious members, are 38% and 18% respectively. 

All these features allow to describe and organize solar 
circle kinematics observations in a simple scenario. 

1. Star formation in the galactic disc occurs in large bursts 




5e+08 1e+09 1.5e+09 2e+09 2.5e+09 3e+09 
Age (yr) 

Fig. 15. Main bursts in the A-F type star sample. Distribution 
of all the stars (full line) and stars involved in streams with 



typical velocity dispersion a stream 
line) 



3.8 km 



( dashed 



separated by quiescent periods (Figure 15). The most 
recent of these bursts started ~ 10 s yr ago (streams in 
Pleiades SCI, IC2391 SCI and Centaurus associations) and 



includes the formation of the Gould Belt. It seems to be the 
start of a new active era after a half-billion year quiescent 
period. Within the look-back time of our sample, there is 
evidence for two long active periods, one around 1 Gyr 
(streams in Hyades SCI, IC2391 SCI, Sirius SCI and New 
SCI) and another around 2 Gyr (the oldest streams). The 
typical duration of burst eras is also of the order of half a 
billion years. It is not clear, given the time resolution of 
this investigation whether during active periods the average 
star formation intensity oscillates on shorter time scales. 
Streams in the 1 Gyr burst reveal preferential ages at 
6 • 10 8 , 8 • 10 s and 10 9 yr favouring this idea of oscillations, 
possibly related to spiral waves. It may also reflect more 
local phenomena. 

2. Under burst conditions, stars mainly form in groups 
reflecting the clumpy structure of the interstellar medium 
sporadically sampled by the star formation process. As 
a consequence, the velocity space is gradually filled by 
successive star formation bursts. Formation puts stars 
preferentially on near circular orbits filling the center of 
the velocity ellipsoid (young streams in Pleiades SCI, 
IC2391SC1 and Centaurus). About 75% of recently formed 
stars belong to streams which internal velocity dispersions 
do not exceed 4 km ■ s _1 . A limited fraction of the initial 
groups are gravitationally bound and form open clusters. 

3. Open clusters sustain a stream of stars with similar velocity 
by an evaporation process due to internal processes or tidal 
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disruption by gravitational potential large scale inhomo- 
geneities (encounters with GMCs). This phenomenon is 
well illustrated by the needlelike density structure found 
around the Hyades OC1. Streams found at scale 1 and 2 
(J? stream ~ 2.3 and 3.8 km ■ s -1 ) seem to be in agreement 
with this scenario despite the fact that age and membership 
accuracy does not always permit to show clear correlation 
between age and velocity. 

4. In this picture the survival of streams as old as 10 9 yr 
is satisfactorily explained as the end of the evaporation 
process of the most concentrated clusters. Observational 
evidences obtained by Wielen, 1971 and Lynga, 1982 



show that the half lifetime of open clusters in the solar 
neighbourhood is between 1-2 10 s yr. A longer survival 
for streams can be explained in terms of resonant trapped 
orbits by the Galactic bar gravitational potential. Indeed, 
all the streams older than 10 9 yr are on the external part of 
the velocity ellipsoid with a high absolute value of their U 
component. They probably have been formed in the inner 
part of the disc where the bar potential is capable to lock 
them into an orbital resonance (Dehnen, 1998). 

5. The typical scale of Eggen's superclusters (a stream ~ 6.3 
km ■ s _1 ) does not seem to correspond to any physical 
entity. For one thing, the picture they form, their fre- 
quency and their divisions at smaller scales are compatible 
with their creation by chance coincidence of physically 
homogeneous smaller scale structures (a s t ream ~ 3.8 
or 2.4 km ■ s _1 ). For the other thing, their internal age 
distributions more or less reflect the overall age distribution 
of the whole sample, with occasionally some preference 
for the typical age of the dominant sub-structure. 



Beyond this phenomenological classification, the 6D analysis 
of this complete sample of nearby A-F stars provides the first 
time dependent picture of the mechanism creating the stellar 
velocity distribution in the disc. 
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